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We derive a geometrical version of the Regge- Wheeler and Zerilli equations, which allows us to 
study gravitational perturbations on an arbitrary spherically symmetric slicing of a Schwarzschild 
black hole. We explain how to obtain the gauge-invariant part of the metric perturbations from the 
amplitudes obeying our generalized Regge- Wheeler and Zerilli equations and vice-versa. We also 
give a general expression for the radiated energy at infinity, and establish the relation between our 
geometrical equations and the Teukolsky formalism. The results presented in this paper are expected 
to be useful for the close-limit approximation to black hole collisions, for the Cauchy perturbative 
matching problem, and for the study of isolated horizons. 



I. INTRODUCTION 



. In the last decade, perturbation theory for black holes has played, in several different ways, a key role in numerical 
' and computational relativity. Already in the seventies it proved to be a very valuable tool to predict gravitational 
waveforms from processes such as a particle falling towards a black hole. Since the early nineties, due to the network of 
04 interferometric gravitational wave detectors in construction, there has been renewed interest in predicting waveforms 
\ for strong sources of gravitational waves such as black hole collisions. In particular, the first predictions using 
' ' I perturbation theory in this new era have been quite striking Some of the applications of perturbation theory in 
' recent years involved computing the evolution for different conformally flat initial data describing black holes in the 
, close limit in order to predict radiated energy and angular momentum [^, to provide both analytical understanding 
' and benchmarking of full numerical results or to quantify the amount of spurious radiation in conformally flat 
T— I initial data Q (see Q for a general review). The usual Regge- Wheeler (RW) - Zerilli and Teukolsky |^ 

formalisms have also been extended to second order ||^, a necessary step in providing first-order perturbations with 
Q their own "error bars" jlO| . Other recent approaches use black hole perturbations to extend the computational domain 
in numerical simulations to the radiative zone via Cauchy-perturbative matching [Tl| ], or concentrate full numerical 
^ resources in the nonlinear regime and let perturbation theory take over in the late stage of black hole collisions . 
bi),. All of the applications just mentioned, though diverse, have a common feature: they are limited to perturbations 
■ of Schwarzschild black holes in Schwarzschild coordinates, and Kerr black holes in Boyer-Lindquist coordinates. The 
, reason for this is that, until very recently, most of the initial data typically used in numerical relativity were for 
maximal slicing, and thus reduced, in the various regimes where perturbation is used (far region, late times, initially 
H close black holes, etc), precisely to the Schwarzschild and Kerr spacetimes in Schwarzschild and Boyer-Lindquist 
■ - - 1 coordinates, respectively. In recent years, however, work was started on Kerr-Schild-type initial data which are 
not maximal. Part of the motivation for introducing this new kind of initial data is to avoid the typical grid stretching 
that maximal slicings produce near the event horizonQ a stretching that eventually causes numerical simulations to 
crash 0. One is then faced with the fact that in order to accommodate these new initial data, either for the close-limit 
approximation or for Cauchy-perturbative matching, a formalism is needed that allows perturbations in more general 
slicings than Schwarzschild and Boyer-Lindquist. 

Another important motivation for having such a formalism in place is to study the recently developed isolated- 
horizon formalism [ pT| in the perturbative regime. For such studies, one needs to be able to analyze a neighborhood 
of the background horizon, which necessitates the use of horizon-penetrating coordinates. 



*sarbach@gravity.phys. psu.edu 
ttiglio@gravity.phys.psu.edu 

^One does not have to move from maximal slicing to get rid of spurious radiation; it is enough to use initial data that is 
conformally Kerr, instead of the more usual conformally flat |Q. 

^ There is some new evidence that these crashes can be avoided by excising the singularity and appropriately choosing the 
shift vector | |l5| . 
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The two most used approaches to black hole perturbations have been the RW-Zerilli and the Teukolsky ones. Each 
of these methods has its own advantages and limitations: The Teukolsky formalism can be used for rotating black 
holes, but one cannot obtain the whole perturbed geometry but, rather, \['4 or (this is enough to compute radiation, 
thoug h)0. The RW - Zerilli technique, on the other hand, provides the whole perturbed metric, but is limited to non 
rotating black holes. 

The Teukolsky equation, in its original formulation, can, in fact, be used to describe perturbations around any 
Petrov type-D background, without relying on a particular choice of coordinates. Work has started very recently on 
the application of this to Kerr-Schild black hole perturbations [|l8| . 

This paper, in turn, develops an appropriate extension of the RW-Zerilli formalism to perturbations of a 
Schwarzschild black holes in arbitrary spherically symmetric coordinates. One can imagine a huge variety of ap- 
plications of such an extension; here we have concentrated on the aspects of the formalism that we need in order to 
proceed with our main motivations. In order to generalize the RW-Zerilli formalism, we start from a perturbation 
formalism introduced by Gerlach and Sengupta pS| ] and derive two master equations which hold in any spherically 
symmetric coordinates of the background, but reduce to the equations obtained by Regge- Wheeler and Zerilli if one 
uses the standard Schwarzschild coordinates. 

Our approach is organized as follows. In section II we present the basic formalism that decouples the field equations 
into the generalized RW and Zerilli ones. The special cases with total angular momentum / = 1 and / = are treated 
carefully. In section III we work out a relation needed for Cauchy-perturbative matching, namely, the one between 
the RW and Zerilli functions and the ADM quantities. In section IV we establish the relation between the present 
formalism and the Teukolsky one, a relation that is desirable not only to compute the radiated energy and make 
contact with |18|, but also from a conceptual point of view. Finally, in section V we comment on the properties of 
the RW and Zerilli equations, and on a numerical code that we have written to solve them. In order to establish 
the contact between the abstract formalism in the body of this paper and more direct applications, we give some 
explicit expressions in appendix A. Finally, in appendix B, we summarize some properties of spin-weighted spherical 
harmonics which are needed in section IV. 



II. THE GENERALIZED RW AND ZERILLI EQUATION 

In what follows, we assume that the background spacetime (M, g) can be represented as a product of M = M/S0(3) 
and with metric 

9 = 9ab dx'^dJ' + 'qab dx^dx^ . (1) 

Here g = dO,^ is the standard metric on S*^, and g and r denote the metric tensor and a positive function, respectively, 
defined on the two-dimensional pseudo-Riemannian orbit space M. In what follows, lower-case Latin indices refer to 
coordinates on {M,g), while capital Latin indices refer to the coordinates d and (p on (S'^,g). Below, we will derive 
perturbation equations which do not depend explicitly on the background metric coefficients. In fact, we will only 
use the background equations which are given by the components of the Einstein tensor 

Gab = - ^ Va Vbr + (2rAr + iV - l) gah , (2) 

Gab = (rAr ~ r^kj cjAB , 
GAb = 0. 

Here, N = g{dr,dr), and k denotes the Gauss curvature of the metric g. A coordinate-invariant definition of the 
ADM mass is given by 

M = ^(1-iV). 



^Actually, one can construct solutions to the linearized vacuum equations from a potential (which is not '^4) that satisfies the 
Teukolsky equation This approach is very interesting but there are some issues that still have to be worked out before 

it can be implemented. For example, how to give initial data to the corresponding potential (in particular, can one obtain 
any linear vacuum perturbation of the Kerr spacetime from some potential?), how to construct gauge invariants and extract 
radiation, etc. 
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We can see by inspection that this M is the mass if Schwarzschild coordinates are used; on the other hand, M is 
defined in terms of scalars on M, and therefore M represents the ADM mass in any coordinates on M. Note that 
in a vacuum spacetime, equation ^ impHes that = r'^{Gab^^'r ~ G'^f^Var) — \7a[r{l — N)] which shows that M is 
constant. 

Since the background is spherically symmetric, it is convenient to expand the perturbed metric in spherical har- 
monics: 

Sgab = HabY, 
SgAb = QfeV^r + KSa, 
SgAB = KgAsY + FVa^bY + 2/cV(^^b), 

where Hab denotes a tensor field, and /if, vector fields, and K , F and k scalar fields on M. Here, Y = F'™ are 
the standard spherical harmonics, and Sa — {*dY)A = e^Vsy and 2\/ (aSb) = ^ aSb + ^ bSa form a basis of 
odd-parity vector fields and symmetric tensor fields, respectively, on (See Appendix D of Ref. pf] for more details 
on spherical tensor harmonics.). We suppress the indices Im and the sum over these indices since the modes belonging 
to different pairs of Im decouple in the perturbation equation. The Y^™ are normalized with respect to the standard 
metric g on 5^, an exception being the cases I = and I — 1: There, we choose the normalization such that — 1, 
and 47r/3. 

In what follows, it will also be convenient to use a coordinate- free notation for differential forms on S and 

S = id* denote the Hodge dual and the co-differential operator, respectively, with respect to g. That is, 

iUadx'^ = iabU'^dx^, 
1 a A T B 

*Vadx = €abv ax , 

where e and e denote the standard volume elements in {M,g) and (5^, 5), respectivel};^. For example, we have 

A further simplification comes from the fact that a spherically symmetric metric is invariant under parity transfor- 
mation X 1-^ —X. As a consequence, the above defined amplitudes decouple into two sets, one set transforming like Y 
(called scalar perturbations or even-parity perturbations) and the other set transforming like S = *dY (called vector 
perturbations or odd-parity perturbations) under parity transformations. In this sense, the amplitudes Hab, Qb, K 
and F have even parity while the amplitudes /if, and k have odd parity. 

A. The odd-parity sector 

We start with the simpler case of the odd-parity sector. The perturbations of g^^ are parameterized in terms of a 
scalar field k and a one- form h = hadx"^ , 

Sgab^O, SgAb ^ hbSA, SgAB ^ ^kV (^aSb), (3) 

where k and ha depend on the coordinates x^ only. Note that for I ~ I, ^(a^b) vanishes and k is not present. For 
I = 0, Sa = and there are no gravitational perturbations. 

1. Coordinate-invariant amplitudes 

A vector field X = X^d^ generating an infinitesimal coordinate transformation with odd parity is determined by 
a function f{x^), where 

X^ = 0, X^ = ^g^^SB. 



*For e we need to provide an orientation in Af; if t and x are timelike and spacelike coordinates, respectively, we choose 
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Using the fact that to linear order, transforms like the Lie derivative of the background metric with respect 
to X, we find the following transformations 

h + r^d( ^ 



k + f. 

Note that one can choose a gauge in which fc = 0. This gauge, which is usually called the RW gauge, is unique. 
For I > 2, one can construct the coordinate-invariant one-form 



For Z = 1, we will see that only the invariant two- form 



enters the perturbation equations. 

In terms of these gauge-invariant quantities, the components of the linearized Einstein tensor are 



l^(inv) 



SGab = 0, 
'Ab I 

SGab 

where the background equations have been used, and where here and in the following, 

A = il-l){l + 2). 



~Jt/j("-)V(^^B), 



(4) 
(5) 



The master equation 



The vacuum perturbations with odd parity are obtained from equation (Q), which yields 



0. 



(6) 



The usual way to derive the RW equation for Z > 2 from equation (||) is to decompose the one- form with respect 

to Schwarzschild coordinates, ft.^™") = /ij™"''(i<-|-/if™"-'(ir, and to use the integrability condition (||) to eliminate /ij™"-*. 
This yields an equation for h^™^ alone, which is then cast into a wave equation for the function <i> ~ {l — 2M/r)h^™^/r 
defined below). This can also be achieved in a coordinate- invariant way as follows: One uses the integrability 
condition = to introduce the scalar potential <& according to /i^™") = id{r^) = eabV°(r$)da;''. Equation 

(O) may then be integrated to yield the following wave equation 



-A- 



$ = 0, 



(7) 



where the two-dimensional Laplacian of a function is A"I> = — d^cZ^ — V"Va$. Here, the free constant in the potential 
$ has been used to set the integration constant to zero. Equation ^ is the coordinate-invariant version of the RW 
equation. Indeed, we have not specified any coordinates on the orbit manifold M. Using the coordinate-independent 
vacuum background equation = r^G"^ — 2{rKr + N — \), equation Q finally assumes the form 



$ = 0, 



(8) 



with 
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Vrw 



l{l + l) - 



6M1 



(,From equation (|^), we also get the following relation 



-id 



(9) 



which enables us to compute $ from the gauge-invariant one-form /i^™^). 

For I — 1 equation (|^) is immediately seen to admit the solution 1/r. Since A = 0, we may also directly integrate 
equation (||). This yields 



id 



6J 



(10) 



where 6J is a constant of integration. At this point, it is important to recall that the one-form h is not coordinate- 
invariant, but transforms according to h ^ h + r'^d{f /r"^). This implies that the solution of the homogeneous part of 
the above equation is pure gauge. A special solution is 



^ _ 2 J idr 

^ 7 17 



(11) 



As explicitly shown in this describes the Kerr metric in Boyer-Lindquist coordinates in first order of the 
rotation parameter a = J /M. By equation ([lO|), J is defined in a coordinate-invariant way. In summary, a general 
I = 1 perturbation is given by 



2J idr 



r^dif/r^ 



(12) 



with / an arbitrary function on the orbit space. 



B. The even-parity sector 

The even-parity perturbations of g^jy are parameterized by a symmetric tensor field Hab, a one- form and two 
scalar fields K and G on the orbit space M, 

Sgab = HabY, 
SgAb = Qb^AY, 

S9AB = KqabY + Gr2 (v bY + U{1 + lYgAsY^ . 

Here, the basis of symmetric tensors in 5gAB is chosen to be orthogonal with respect to the inner product induced 
by g. Furthermore, one has V a'^ bY + + ^)gABY = for Z = 0, 1; hence the amplitude G is not present in those 
cases. For I = 0, the amplitude Qb is also absent. 



1. Coordinate-invariant amplitudes 

An infinitesimal coordinate transformation with even parity is generated by a vector field X with 

X'' = i^^Y, X^ = fg'^^^BY, 

where and / are a vector field and a function, respectively, on M . With respect to this, the metric perturbations 
transform according to 

Hab ^ Hab + £,a\b + 6|o, 

Qb^Qb + ib + r''f\b, (13) 
G^G + 2f. 
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Here and in the following, = '^a^b denotes the covariant derivative with respect to the orbit metric g and 
Va = r\a/r. 

For I > 2, one can construct the following set of coordinate-invariant amplitudes: 

H^r^ ^Hab-{Pa\b+Pb\a), (14) 
J^{^nv) ^ ^ _ 2t,V + + (15) 

where Pa — Qa — \r^G\a- The (generalized) RW gauge is defined by choosing and / such that Qi, and G vanish. 

We see that in this gauge, which is also unique, and ift"'') coincide with Hab and K. 

For 1 = 1, there is no such simple choice of coordinate- invariant amplitudes, since G is not present in this case. 
Nevertheless, we can always chose the gauge such that Qb vanishes. One then remains with Hab and which are 
subject to the residual coordinate transformations as in (^3|) with ^f, -|- r^/|j = 0. 

For I = 0, Qb and G are absent anyway and one can arrange the gauge such that K = 0. 

In summary, it is sufficient to derive the linearized Einstein equations for the perturbed metric 

Sgab = HabY, SgAb = 0, SgAB = KqabY, 
where for I >2, Hab and K can be replaced by their coordinate- invariant counterparts defined in (0,f5|). 



2. The master equation 

The long but straightforward computation of the linearized Einstein tensor is given in [20| . The equations' structure 
becomes much more transparent if one first splits the two-tensor Hab into its trace and traceless part and then 
introduces the one-form 

C = Habr\''dx\ (16) 

where Hab denotes the traceless part of Hab- A similar split is performed for the components of the Einstein tensor. 
As a result, the relevant components of the Einstein tensor define two scalars S and T and two one-forms U and V 
according to 

SGab = S (vAysY + h{l + lygAsY^ , 

r'SGab^TY, 

SGAbdx' = ^U\7aY, 

SGabT^'^dx'' = VY. 

The vacuum field equations are then expressed in terms of the one-form C and the two scalars H = g°'^Hab and K. 
The simplest equation, which is present only for I > 2, gives 

= -25* H, 

hence Hab is traceless (For I = 0, 1, we can make use of the residual gauge freedom in order to impose H = 0. Residual 
coordinate transformations are then of the form (^ with = and = -r'^f\a for I = 1.). Using H = 0, the 
remaining equations reduce to 

= T = - Jtc - ^g(C, dr) +Kk+ -gidK, dr) - ^K, (17) 
= f/ = -^ {d1C)dr + {idG)idr - dK, (18) 



= V= [d^G)— + -d~g{C, dr) + ^^^^r^C 



2r2 

, , , Ar 

2 , . , 2r 



^KKdr- dg{dK, dr) + ( Ar - \ dK, (19) 
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where equation is void for / = 0. We recah that for Z > 2, we should replace Hab by H^^^^ and K by JT^™") 
in the above equations in order to give them a coordinate-invariant meaning. For ^ > 1, we compute the component 
of U parallel to *dr: 

= giU, idr) = idC - g{dK, idr) = Sd [C - rdK] . 

This motivates us to replace the one-form C with the one-form 

Z = C - rdK. (20) 

In terms of Z and K Einstein's equations become 

= g{U, idr) = idZ, (21) 
= g{U,dr) = -d'fZ + rAK, (22) 

= T = -dtz - ^g(Z, dr) -~ Kk ~ irK, (23) 

= {2V - Tdr) ^ d [2rg{Z, dr)] + l{l + l)Z + r(ao + \)dK + XKdr, (24) 

where we have defined 

ao = 2rAr + 1 - TV. 

Using the background equation = r^G"^^ = 2(rAr + N — 1) and = 1 — 2M/r, one finds ap = 6M/r. 

In view of equation (pl|), we may introduce the scalar field ( according to Z — d(. Equation (|2^) may then be 
integrated to yield 

2rg{dC, dr) + l{l + !)( + r(ao + X)K = 0. (25) 

It is now clear how to obtain a single, second order differential equation for C: First, we eliminate AK from equations 
(H) and (|2|). This gives 

- AC - -g{dC, dr) --K = 0. (26) 
r r 

Next, this equation is used to eliminate K in (|25|). Hence, 

- (ao + A)AC - dr) + ^^^^C = 0. (27) 

(Note that for / = 1, this equation is equivalent to equation (^) and thus is also valid in that case.) Finally, we define 
the new scalar function by 

C - (ao + A)^', 

in order to remove the first order derivatives. This yields the Zerilli equation Q], 



-A + Vz 



= 0, (28) 



where 



_ A^r^ [(A + 2)r + 6 Af )] + 36AP ( Ar + 2M) 
^~ (Ar + 6Af)2r3 ' ^ ' 

Before we discuss the special cases I — and ^ = 1, we make two remarks: First, for I > 2, the scalar field 'J can 
be obtained from the Zerilli one- form Z using equation ( p7| ) . The second point is that it is also possible to obtain the 
RW equation for the scalar <I>e = r'^d} Z. In fact, Chandrasekhar (see, e.g. |27|) has shown that the equations of RW 
and Zerilli for a Schwarzschild background are equivalent in the frequency domain. However, in the time domain, we 
were not able to express Z in terms of $6 and its derivatives alone. For this reason, we will use the Zerilli equation 
in the even-parity sector and not the RW equation. 

For I = lj equation ( p7| ) reduces to 
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- AC - ^g(dC, dr) = ^d^r'dC) = 0. (30) 

However we recall that for I — I, Hat , and hence also C are not defined in a coordinate-invariant manner. Under 
the residual gauge-freedom (|l^) with = —r^fia a-nd 

dt(r2d/) = -d" = 0, (31) 

we find that Z transforms according to 

Z 1-^ Z + 6M df, and hence C C + QMf. 

Since / is an arbitrary solution of equation (^TJ), and since the equations (|30| ) and (31) are equivalent, it is clear 
that every solution of ( pOl ) corresponds to pure gauge. In particular, we can choose the residual gauge in order for ( 
to vanish. In this gauge K vanishes as well, as a consequence of equation (^5|). The even-parity sector is therefore 
empty for / = 1. 

For I = 0, one can choose the gauge such that both H and K vanish. Then equations ( p7|) and (|l^) yield 

Jt(rC) = 0, g{rC,dr) = 26M = const., 

which has the general solution 

C = '^^dr + idh. (32) 

Here, 6M is a constant describing the variation of the ADM mass, and /i is a function that only depends on r. 
Comparing this with a residual gauge, which is generated by — eatk^'' for a function k of r, we get 

C ^ C ~ i{Nk"{r)dr), 

showing that the function h{r) above corresponds to pure gauge. This can also be seen in a gauge-invariant way: 
Recall that for any spherically symmetric metric of the form (|l|) we defined the mass parameter M through 1 — 2il//r = 
N — g°'^r\ar\b ■ Using the fact that (for Y=l) Scjab = Hab and (5(r^) = r^K, we obtain 

25M = rr^^'r^^'Hab - rr^''irK)\b + MK. (33) 

It can be checked that the RHS is indeed a gauge-invariant combination. On the other hand, for K — equation 
(^ yields 25M = g{rC,dr), as above. 



C. Summary 

In both the odd- and the even-parity sector, perturbations on any spherically symmetric vacuum background are 
described by a wave equation of the form 

(-A + y)u = 0, (34) 

where A is the Laplacian with respect to the orbit metric g and where the potential V depends on the ADM mass 
M, r and the angular momentum number I only. 

For I = and I — 1, there are no dynamical modes. The only physical solutions in those cases are stationary, 
describing variation of the mass and angular momentum. The gauge-invariant part of the metric can be obtained 
from u and vice-versa. These relations are going to be made more precise in the next section. 

Finally, we would like to mention that our gauge-invariant perturbation formalism has also been generalized to the 
case where matter fields are coupled to the metric [ pOpl| ]. In the case of Einstein-Maxwell, we were able to generalize 
the equations obtained by Moncrief ||2^. However, as we have argued in a recent Letter the perturbation 

formalism presented here fails to yield a wave equation of the form (p4h with a symmetric potential, V — V'^ , when 
non-Abelian fields are coupled to the metric. 
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III. RELATION TO THE ADM QUANTITIES 



As mentioned in the introduction, one of the motivations for the present work is Cauchy-perturbative matching. 
This amounts to matching numericaUy, at each time step, the variables used in a nonhnear code with the ones used in 
the perturbative regime (in our case the RW and ZeriUi functions). The matching takes place at a timelike boundary. 
For this purpose we explicitly show the relation between the RW and ZeriUi gauge-invariant potentials and the ADM 
quantities, namely, the three-metric and the extrinsic curvature. This does not restrict the formulation of Einstein's 
equations to be used in the nonlinear regime, since for a formulation other than the standard ADM (e.g. conformal 
ADM, or a hyperbolic formulation) the relevant quantities can be obtained from the three-metric and the extrinsic 
curvature, and vice-versa. 

So our aim is to make explicit the relationship between the scalar fields $ and 5" satisfying the RW and ZeriUi 
equations (||) and ( p8| ) and the components of the linearized 3-metric and extrinsic curvature. We will show in this 
section that - modulo gauge transformations - there is a one-to-one correspondence between 6gij, SKij and the scalar 
amplitudes $ = 9t$, 5*, 5*. Furthermore, this correspondence involves no time-derivatives. For example, it is 
possible to express 4' in terms of purely spatial quantities, i.e. Sgij, SKij and their spatial derivatives only. 

We assume that the full metric, satisfying the nonlinear field equations, has the ADM form 

gr(^) = -a{fifdt^ + g^jifi) (dx' + P\y:)dt) {dx^ + l3^fi)dt) , 

where u is a variational parameter, such that for /i = 0, the metric is spherically symmetric. With respect to the 2+2 
split (|l|), the orbit metric g takes the form 

g = -a^ dt'^ {dx + (3dtf , 

where x is any radial coordinate, a and f3 = are the background lapse and shift, respectively, and 7^ = gxx- The 
components of the extrinsic curvature are 

2aKxx = 27 (9o7 - iP') , 
2aKxA = 0, 
2aKAB = 2rdorgAB, 

where a prime denotes differentiation with respect to x, and where we have also introduced the normal derivative 
do = dt- Pdx- 

The components of the linearized metric have the form 

6gtt = ~2a6a - p^Sg^x + W^Px , 

Sgtj = S/3j, (35) 

where i = x, A. Note that we use perturbations of the coshift rather than the shift vector. This fact will turn out 
to be important when we try to express the ADM quantities in terms of the RW and ZeriUi scalars. Similarly, the 
components of the linearized extrinsic curvature are given by 

2aSKxx = dtSgxx - 2KxxSa + ( ^] - 2j 



j2 J . y ^ 



6/3 

2aSKxA = dtSgxA - 2(3—SgxA + A^gxx - '^aS(3x - ( 



2a5KAB = doSgAB - 2KABSa + 2/3V^ASgB)x ~ 2W(aSPb) (36) 



rr' 



'^—^9AB {S(3x - PSgx. 
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A. The odd-parity sector 



In the odd-parity sector with I > 2, the only non- vanishing perturbations can be parameterized according to 
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5f^A = bSA, 

SgxA = hiSA, 5gAB ^ 'ih2'S/ (aSb), 
^Kxa^t^iSa, 3Kab = "27^2^ (aSb)- 



1. The potentials in term of the ADM quantities (I > 2) 

We want to express $ and $ in terms of the quantities b, /i2, tti, tt2 and their spatial derivatives. 
First, we observe that ht — h, = hi, and k — h2 where ht, and k are the ampHtudes introduced in 
Therefore, we obtain 



hr^^ ^b-h2 + 2-h2 (37) 
r 

= /ii _ h'^ + 2^-h2 . (38) 

Next, the equations ( ^6| ) yield the relations 

r' / b\' 

2a7ri = hi- 2(3-hi - rM ^ , (39) 



r 

2an2^dQh2+Phi~b. (40) 

Using the last of these relations, eq. (^), to re-express, in eq. (^), time derivatives of /i2 in terms of spatial 
quantities, the components of the gauge-invariant one-form take the form 

4™")= -2^712 2^/12, (41) 



r 



4"") =h,-r'{^] , (42) 

= /if ™) - Next, on( 

(|39|), to find 



where /ig*""^ = /if"^^ — (3hx™^ ■ Next, one uses equation (p|), trading time derivatives for spatial ones with the aid of 



$ = — ( 2a^i -2^/ii ) , (43) 



Aa7 



r 



which is one of the formulae we were looking for. In order to obtain the time derivative of one uses the definition 
of $, i.e. /i^™") = ;rf(r$). This yields 

4-) = _^9,(,$), (44) 
7 

h^^''^ = -^do{r<^). (45) 
a 

which one can solve for <i>. Using eqs.( |4^j4l| , ^ ), the result is 

• 1/ 2fdor\ 2^13 2r ra ( h2\ 2^(3f 

$ = — -a H r hi H 7r2 - — tti H — — /12 (46) 

7r \ Aa J r A7 1 / ^ 



2. The ADM quantities from the potentials (I > 2) 

On the other hand, given $ and $, we obtain Sgij and 6Kij in the following way: First, we compute /iq""'' and 
^(inu) equation ( ^4p^ ). Then, using the above equations, it is straightforward to express b, hi, tti and tt2 in 

terms of h^Q™\ hx^^\ $ and k, where k parameterizes the gauge freedom: 
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b^ht^^+(3h(:"^^+r'dt(^^y (47) 
h, = /i^"") + r^d, (\) , (48) 



h2 = k, (49) 

2a^i = + 2— /ii , (50) 

r r 

2a7T2 = -hl"'^'> +2^k. (51) 



The special case 1 = 1 



For / = 1, the amplitudes /12 and 1:2 are absent. According to the analysis in the last section, the only physical 
solution is the Kerr mode. Using equation (|lO|), one finds that the rotation parameter (the only gauge invariant for 
Z = 1 ) can be extracted from the ADM quantities according to 

6J = — (2aTTi-2—hi\ . (52) 



aj \ r 

On the other hand, using (|ll|), one finds 

2J //37^ a ,\ , 2^ // 



, 2J7 ^ 2a ff\ 

ri\ a X"^ J 

where / parameterizes the gauge freedom, and where = 1 — 2M/r = —{dorY'/cP' + The amplitude tti then 
follows from (p3). 



B. The even-parity sector 

Here the perturbations are 

5a ^aY, S(3, = biY, 5f3A = b2VAY, 
Sgxx ^ hY, SgxA = Q^aY, 

6gAB = KqabY + Gr^ (vaVbY + + IjgABY ] , 



SKxx = TTh Y, SKxA = ng\7AY, 

7t,Y -4- ^„ 



SKab = TTK9ABY + TTG rM VaVbY + h{l + l)gABY 



1. The potential in terms of the ADM quantities (l>2) 



IYtotcl equations (35) one finds Hu = —2a a — f3 h + 2/35i, Hfj. = 61, Hj.x. = h, Qt = b2, Qx = Q, while K and 

The expressions for the linearized curvature tensor, equation 



G agree with their definitions in the previous section. 
(§), yield 
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2a = /. + ;37M ' - 27 (^) ' - ^ (90 - 7/?') , 
2aTrg = q-2l3'-q + f3h-bi-r^ (^^^ , (54) 
2aTrK = \do{r^K) + ^(fih- h) - ^-^^^ {(3q - &2) - —dor a, 

1 „ . 2 2 



At first sight it is not clear how the ZcriUi one-form Z, defined in ( |20| ) can be expressed in terms of spatial 
amplitudes only, since from the definition of i?^™"^ one sees that second time derivatives of metric components can 
appear. However, it turns out that only the two-form ujab = Pb\a — Pa\b^ which contains no second time derivatives of 

h, q, K and G, appears in the Zerilli one-form. Using the fact that h'^™^ is traceless as a consequence of the field 
equations, one obtains 



Now, using ([5^), it is easy to find 

Po = -ar^TTG + r{dor)G, (55) 
P. = q- \r^G', (56) 

c^o, = 2a{nq -f r/TTc) - ^{dor){r^G)' - (/3/i - 61). 
Using the background equation v^a = Mr~^gab — f'~'^r\ar\b, one eventually obtains 

Zo = -2ar (ttk + h{l + 1)^g) + 2(9or) (k + h{l + 1)g\ (57) 



dor \ 2 ( M 



2a-Kq - 2—q - - 1 - 3— po 



1,,, ,,^V ^0?- c^dar \ 2 / 



Z, = -^h-r[K+ + l)Gj + ^ [2a^q ~ 2^qj - ^ - ) ■ (^8) 
The Zerilli scalar C (and ^) can now be obtained from its definition, eq. (p5|) , with 

/^(--) =if + iz(; + l)G--rl%, 
2 r 

r^''pb^ -{dor)po/a^ + r'p^/j^ , 

and the one- form Z given by equations (|5^,^). On the other hand, the latter equations also give us from = 
Zq + /3Zx ■ Note that - as in the odd-parity sector - the scalars ( and C do not depend on the perturbed lapse nor on 
the perturbed shift. Finally, we see that for a Schwarzschild slicing where dor = 0, $ and ( are linear combinations of 
the extrinsic curvature components only. These combinations precisely agree with the ones obtained in a perturbative 
approach on a static background in terms of curvature-based quantities . 

2. The ADM quantities from the potentials (I > 2) 

If C and C are known, equation ( p5|) tells us how to obtain X*^*"^-' and K^'-™') . Next, the traceless part of i?^™"'' is 
obtained from this and the definition of the Zerilli one- form Z . Finally, one has 
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2aa = -F^r' - 2^,^ + 2/?(pt|, + p,,*) - 2/3^p,|„ (59) 
h = HiT^ + 2p,|„ 

^1 = -f^tr"^ +Pa;|t, (60) 

b2=Pt + \r^G, (61) 
r ' 2 

Here, £ind G parameterize the gauge freedom. The amphtudes tt^j, ... ttq are obtained from this and (|5^). 

5. r/ie special case 1 = 1 

For Z = 1, according to the analysis of the last section, we have only gauge modes, and G and ttq are absent. 
Therefore, the ADM-based amplitudes are obtained from the same equations as above, but where H^^^^^ and K^"™^ 
are set to zero, and Pa and G are replaced by and 2/, respectively, where and / parameterize the gauge 
transformation that brings us from the RW gauge to the actual gauge that one wants to use. 

^. The special case / = 

For / = 0, 62, 9 and tt^ are also absent. Using equation ( p3| ) and the relations ([5^), the perturbed mass parameter 
is found to be 

5M = '^{d,r).K + \{^)'h+\{r- M)K - ^.{r^K)'. 

In order to obtain the perturbed three-metric and extrinsic curvature in terms of 6M, one uses equation (^) which 
gives 

4(5Af / N \ 

= - y5a& j + gauge, 

and the ADM quantities are obtained in a similar way to above. 

C. Gauge fixing vs choices of lapse and shift 

We have shown above how to construct the three metric and extrinsic curvature from the potentials (and vice- versa), 
up to gauge freedom. In numerical simulations, however, usually one does not fix the gauge but rather chooses lapse 
and shift, perhaps as prescribed functions of spacetime ("exact lapse" or "exact shift") or as dynamical quantities 
coupled to the the three metric and/or extrinsic curvature ("live gauges"). In general this does not fix the gauge 
completely, which means that we have to relate the gauge freedom to the choice of lapse and shift. The properties 
of such relations depend on the details of how the lapse and shift are chosen, and it is therefore not possible to give 
a general discussion. These equations, for example, might be elliptic if some kind of minimal distortion is imposed, 
hyperbolic as in the case we discuss below, or of some other (perhaps unknown) type. 

Here we will concentrate on a specific simple prescription, but it should be clear that other cases can be treated 
similarly. The case we are going to discuss is exact-coshift, exact-lapse; that is, the lapse and shift covector are 
arbitrary but a priori given functions on the orbit space. 

We start with the odd-parity sector. The perturbed lapse is zero, and, for I > 2, the perturbed coshift is given by 
the right hand side of equation (|4^). The function k parameterizes the gauge freedom, and it is thus given by the 
equation 

dtk = ^ {b - hi""" - ph^") (62) 
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Since 6 is a given function, this equation can be solved for k, provided we supply initial data. Given any three-metric 
and extrinsic curvature at to, the initial data for k is given by k{to,x) = h2{to,x). 

The treatment for ^ = 1 is similar. Now the gauge freedom is controlled by /, which can be related to the coshift 
by equation (p3), rewritten as 



■ 2J //?7 a , 

b — dor + -r 

ri\ \ a 7 



(63) 



In the even-parity case with I > 2 the gauge functions pt , Px and G are related to the lapse and shift by evolution 
equations which are straightforwardly obtained from equations ( |59| , |60| , |6l| ) . These evolution equations form a 3 x 3 
coupled system, first order in space and in time, 

G = + 1.0. (64) 
Px = -p't + l-o. (65) 
Pt = -f3^p', + 1.0. (66) 



where l.o. stands for lower order terms. Initial data for the system ( p4y6qj66| ) is given by the three-metric and the 
extrinsic curvature at some time to and the formulas ( p5| , p6[ ) for pt and Px- It is easy to see that equations ( |64| , |65|j6^ ) 
constitute a weakly hyperbolic (see, eg, |^J]) system if /? = and a stricly hyperbolic system otherwise. That is, if 
these equations are written as Ut — Aux+ l.o., with u = {G,Px,Pt)'^ , the matrix A has three different real eigenvalues 
if /3 7^ 0; and a single degenerate real eigenvalue (zero) with only two independent eigenvectors if /3 = 0. The structure 
of the equations for / = 1 is the same, replacing pa, G by ^a, 2/, respectively. Finally, for / = the situation is similar 
but simpler: G does not appear, and the principal part of the evolution equations for two gauge quantities Px and pt 
is also given by (|65|,|6^). As before, these equations are weakly hyperbolic if the background shift is zero, and strictly 
hyperbolic otherwise. If we use densitized lapse, as is usually done in hyperbohc formulations (see, e.g., |Q), the 
above system of equations is strongly hyperbolic even in the case where /3 = 0. In contrast to this, the system is ill 
posed if we use exact shift instead of coshift. For I — 0, this fact has already been noted in Efl]. 



IV. RELATION TO THE TEUKOLSKY FORMALISM 



In order to compare our perturbation equations with the Teukolsky equation for a non-rotating background, we 
introduce a NP null tetrad that is adapted to the spherically symmetric metric (^, i.e. 

I — ladx"", n — Uadx"" , m = niAdx^, 

where I and n form a null dyad of g, 

gab — —lanb ~ h^a , 

and m is a complex one-form such that 

r^gAB = mAffiB + mBfriA ■ 
Here and in the following, a bar denotes complex conjugation. Note that 

Eab = laTT-b ~ , r^^^AB = « {mAfriB - mBfriA) ■ 
The only non-vanishing NP coefhcients are 

1 ^ 1a 

p = Dr, /i ~ — Ar, 

r r 

7 = ^dn{l,n) = -i^VVa^fc, 

a = —(3 = —a, 
r 
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where a = —^dm{m,rn) = ^rh rh^V a^tt-b is a NP coefficient with respect to the dyad defined by m = im. Here, 
D = /°Va and A = n°Vo ■ We also introduce, for later use, the angular derivative operator S = iji^S/a ■ From the 
NP vacuum equations (see, e.g. [^), it then follows that all Weyl scalars but '^2 vanish, ^0 = ^1 = 'I's = ^'4 = 0. 
In terms of the invariant definition of M given in section 2, 5*2 can be expressed as 



^2 = -- 



M 



In particular, the metric is of type D with repeated principal null vectors aligned with 1°" and n°. 

In what follows, we study the decoupled equation derived by Teukolsky [p| governing linear fluctuations of ^'4 
on any spherically symmetric vacuum background. To linear order, ^'o and '^4 are invariant with respect to both 
infinitesimal coordinate transformations and null tetrad rotations. The reason why we focus on ^'4 and not ^'q is that 
we want to study outgoing radiation at null infinity, which is described by '^^ (see [^). However, by interchanging I 
with n and m with to in what follows, one easily obtains the corresponding results for ^Pq, describing ingoing radiation 
at the event horizon. 

With respect to the chosen null tetrad, the pulsation operator acting on the linearized field 5*4^ •* splits into the sum 
of an orbital and an angular operator. 



A 



0, 



(67) 



where 



= (A + 27 + 5//) (D + Ae-p)- S^-a, 



A= - l5~ 2a] [S + ia 



(68) 



Next, we compute the perturbed Weyl scalar ^'4^'': With respect to the background metric (IT]), one obtains 



^4 - ^AbCd 



Performing the multipole decomposition as described in section ||, we obtain 



(1) 



{inv) 



m^Th^VASB 



in the odd-parity sector and 



4- 



(1) 



rTi^TO^VAVpF 



in the even-parity sector. Using the definition of the derivative operator 5 and the NP coefficient d, one can check 
that in both parity sectors, the angular part is proportional to the spin-weighted spherical harmonics Y"!™ defined in 
appendix H. Explicitly, we have 



Y_ 



Irn 



1 



5-2a]5Y 



rlra 



where Cf — {I — 1)/(Z + !)(/ + 2)/4. It remains to express /i^*"^ and i?^™" in terms of the scalar fields $ and C,: 
Using the definitions of the RW and Zcrilli potentials $ and C, as well as equations ( pO| ) and (|2^), we obtain 



1 



ao + A 

where we recall that oq = 6M/r. Eventually, we get 

1 



(1) 



E 



r(ao + A) 



(A -I- 27 2/Lt) A(ao + A)^' 



Im 



(A -t- 27 -I- 2/i) A$ 



Im 



2 1 



(69) 
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which gives '^'^'^ in terms of the RW and ZeriUi potentials and Vf/m = Cim/(ao + ^) introduced in section 2. Here, 
we have re-introduced the indices I and m. In order for the metric perturbation to be real, we must have — ^i-m 
and = ^i-m (a bar denoting complex conjugation). Equation (M) (or its Fourier transform in time) is some kind 
of generalization of the Chandrasekar transformation (see , also [|30[ ) . 

The corresponding expression for follows after the replacements A ^ D, j — e, /i ~p and Fi™ t-^ Y.j"^ in 
the equation above. 

Provided that $jm and "^im satisfy the RW and Zcrilli equations (|^) and (^8|), respectively, and using the vacuum 
NP equations and the commutation relations 

iD + 2{s + l)e + qp) {A + 2s-f + p p) - {A + 2{s-l)-f + p p) {D + 2se + q p) = + 2(.s + p + g)*2 , 

where s, p and q are arbitrary real numbers, one can show that indeed, ^'4^'' satisfies the Tcukolsky equation (|67|). 
Finally, the total radiated energy per unit time can be obtained from 



dE r 
—— — lim — , 

du r^oa 47r J g2 



'^^{u, r, rt)du 



dn , (70) 



where asymptotically flat coordinates and an asymptotically flat NP tetrad are chosen, and where u = t — r. In 
our case ^'4 = on the background, so the radiated energy depends only quadratically on . Since the fields 
^im and are scalars with respect to the background metric g, we can evaluate ( |70| ) using any asymptotically 
flat coordinates on the background. Using the fact that at infinity, A = \{dt — dr) + 0(r^^), a,^,p = 0(r^^), and 
imposing the outgoing wave condition + ^J^^^ = 0, ^;„j + = at infinity, one arrives at 



l>2 m=-l ^ ' 

(In the derivation, we have also used the orthogonality of the Fi™ and 4';„j = '^i-m , ^im — ^i^-m ■) As a consistency 
check, it is useful to note that this coincides with the usual well known result for Schwarzschild black holes in 
Schwarzschild coordinates |^ Eq.([7ll), however, holds for any coordinates. 

V. THE RW AND ZERILLI EQUATIONS AND NUMERICAL ISSUES 

The RW and ZeriUi equations are wave equations with exactly the same differential operator (the Laplacian on 
the orbit space); they differ only in the corresponding potentials, cf. eq. (^4|). This simplifies their analysis, since 
properties such as well-posedness do not depend on lower order terms (as the potentials are). We now discuss certain 
properties of the RW and ZeriUi equations. In particular, we wish to note the fact that they are perfectly well defined 
as long as the background is regular (both in the coordinate and curvature sense) . Provided the latter holds there is 
no pathology in the equations (or the solutions) at, for example, the event horizon. 

We start writing these equations explicitly by introducing coordinates. We then express the whole metric as 
gt^tai _ gg^^ ^ with the background metric given by 

g = (-a^ + -f^P^)dt^ + 2j^l3dtdx + -f^dx^ + r^{d§^ + sin^ (72) 
The RW and Zcrilli equations are 

Z = ciZ' + C2Z" + C3Z + aZ' - o?YZ (73) 
where Z denotes either the RW or ZeriUi functions. The coefficients c, are 



^Taking into account, of course, that different normalizations are used in the literature when defining the RW and Zerilli 
potentials, see, e.g. |pl[. 
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Cl = 2/3 



C2 = 



C3 



(70; — ^Pa' + q;/37' — wy + jaP') 



7a 



7-^0! 

and the corresponding potentials are 



Vrw — 



r 



Vz = 



A2r2[(A + 2)r + 6M)] + 361\yP(Xr + 2M) 
(Ar + 6Af)2r3 



To make the hyperbohc character of these equations manifest, we introduce new variables y :— Z' , w :— Z, and 
write it as a first-order system for the "vector" u = {w,y,Z)'^, i.e., u = Au' + Bu, where in this case the principal 
part is 



A 




and has eigenvalues and eigenvectors 



Aq = , with (=0 
1 



A+ = 



[0,0,1] 
ci±(c?+4c2)i/^ 



with e± — [A±, 1, 0] 



In our case, cf + 4c2 = 4,a'^j~'^ and, so, the eigenvectors of A are independent provided the background metric 
is locally well defined. Thus, the system is strongly hyperbolic, which is enough to prove well-posedness for the 
initial-boundary value problem if one gives boundary data for the characteristic modes that enter the domain [ p4[ . 
For the close-limit evolution of black holes in horizon-penetrating coordinates, for example, one would put the inner 
boundary inside the black hole, check that the characteristic modes are indeed leaving the computational domain 
(i.e. that the eigenvalues of A are positive), and thus not put boundary conditions there ("excision"). At the outer 
boundary one would typically put zero boundary conditions for the ingoing modes. 

One of the additional advantages of having a hyperbolic equation is that one can write codes that can a priori 
be shown to be convergent |^^. We have indeed written two such codes for the RW and Zerilli equations with an 
arbitrary background. One of them uses fourth-order centered differences in space and fourth order Runge-Kutta in 
time. It uses extrapolation at the inner boundary (assumed to be inside the black hole), and gives zero boundary 
data to the characteristic mode that enters the computational domain at the outer boundary. The other code is 
second-order; it also uses Runge-Kutta for time integration and centered differencing in space, but now needs some 
dissipation (one can prove that this scheme is unstable without dissipation, see js^). In future work we will present 
numerical details of these codes applied to the close-limit collision of superposed Kerr-Schild and Painleve-GuUstrand 
black holes. 
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APPENDIX A: PERTURBED FOUR-METRIC IN TERM OF THE POTENTIALS 



Here we will give explicitly some of the expressions used in the body of the paper, 
coordinate system for the background metric, as in eq. ([72|). 



That is, we choose a general 
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1. Odd-parity sector 



a. Four metric 



The perturbation for the four metric with Z > 2 is given by (Y^ = d^jX , etc.) 



5gm 



a \ J r 



sini^ 
sin dY{) 



2k 
sin^ i9 



' cos dY^ + sin -dY^^] 



k [cos + sin ^ i^Y^^ — sin dY^^] 



Sg-at 



"fa V / r 



sini? 



5^00 = 2k [cos'dY^ - sini9y^0] 



7a V / 



rk — 2kr 



sin 'dYi) 



It is straightforward to compute the Hnearized Ricci or Einstein tensor for the above perturbed metric and see that 
they are indeed annihilated if the master equation (j?^) holds. 

The I = 1 components of the metric, on the other hand, are given by 



f'r - 2fr' 2Jj{-r + fir') 



Ya 



5gx4 



r (2Af ~r)a 

f'r - 2fr' 2 J7(-r + /3r') 



sin 1? 
sin -dY^ 



5gt4. = 



r (2M - r)a 

fr - 2fr 2 J {j^f3{(3r' - f) - a'^r') 
r ^a{2M — r) 

fr - 2fr ^ 2 J {-f^P{f3r' - r) - aV) 



Y^ 
siniS 



sin dY^ 



"fa(2M - r) 

and it is also straightforward to check that this linearized metric satisfies the linearized vacuum equations. 



6. Perturbed ADM quantities 

The three metric can be obtained straightforwardly from the spatial components of the four metric explicitly given 
above and, similarly, the coshift can be obtained from 6 Pi = Sgu and the above expressions for Sgu- On the other 
hand, the nontrivial components of the perturbed extrinsic curvature can be computed directly from the four-metric 
above explicitly, or from the results in the body of the paper. In either case, the results for / > 2, are 



SKxS 



Y, 



4> 



r'^cP- sini9 



a[r - Pr'){k'r - 2kr') + -ir'^{~r + - /3$')+ 



y {-2{r'Yp'^ + a^{l{l + 1) - 2) - 2f'^ + ifirr') $ 
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¥,1 sin 1 



[air - (3r')(k'r - 2kr')+ 



7r2(-f + (3r'){^ - /?$') + ^ {-2{r'ff}^ + a^{l{l + 1) - 2) - 2f^ + 4f}rr') $ 



and for I = 1, 



+ r'$) + — i 

7 ar 

1 [1/ ^/ /^N 2k{r-(3r') 
- -(r$'+r'$ + — 

2 [7 ar 



cos 'd Yj 



+ 



i90 



sin^ 1? sin 

+ cos 'dY^ — sin 



-(r$' + r'$) H — ^ 

7 ar 



\sin^? 
{Y^ cos ■!? — sin dY^^) 



6K. 



{r - pr'){fr - 2fr') 



(472M/366' - 272Ma2 - 272Mr2 - 2/3%^M(6')^+ 



r7' 



r7a2(r - 2M)2 



sint? 



{f-(3r'){f'r-2fr') 



na'^{r - 2M)2 

'a^ - AP^H'bb + 2p'^j^{b'fb - Aa'^{r'fr + &a'^{r'fM + 2(7^)%^) 



(4-f^Mj3bb' - 2-f'^Ma'^ - 2j^Mr^ - 20^-i'^M{b'f 

Y^ sint? 



2. Even- parity sector 

The expressions in this sector are also straighforward to obtain, but the final expression are too long to be written 
down here. For this reason we wil only present the simplest explicit example; the perturbed four metric, in the RW 
gauge, for Z > 2 perturbations of the Painleve-GuUstrand spacetime. The background metric is, thus, given by 



(2MV'^ 
9xx = ^ , 9xt = ( — I 9^^ = r 



= r sm tf Qtt = — 1 H (where r = x) 

r 



while the perturbation is 

f 2MY'^ f rUl + l) -2r + 2M 



^ _ 2^^" _ (^_ r{r + 3M)l{l + l)-2r^^6Mir-My , ^, 



Y 



SQx 



^\ r ) \6M -2r + rl{l + 1) J ^ ~' " \~ r{6M - 2r + rl{l + 1)) 

+ (3r3/5 + r^l^ - ll'^r^ - l^r^ + Alr^ - ISMr^f + UMr'^fi + QMrH'^ - 2AMr^{l - 1) + iQM^rl{l + 1)+ 

72M2(M-r)) 
{QM -2r + rl{l + l)fr'^ 

\-2r^' - ,rir-3M)lil + l)-2r^ + 6Mir-M) . / _ . 

r3/2(6M-2r + rZ(/ + l)) V ^ ' )^ 

(2M)i/2 {2,rH^ + r^f - ifr^ - l^r^ + 4lr^ - ISMrH'^ + l2MrH^ + QMr^l^ - 2AMr^{l - 1)+ 
36M\l{l + 1) + 72M2(M - r)) 



(6M-2r + W(/ + l)) r5/2 



Y 



6g.. 



-{2rf/^M^/H - 2{r - 2M)r*' 



-12M(r - 2M) + QMrl{l + 1) - r^P + 2rH^ + rH"^ - 2rH" 
6M-2r + W(/ + l) , 



Y 
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/„^V2,.l/2i ^/ X/ f -12M (r ~ 2M) + 6Mrl(l + I) - r^P + 2r^P + T^l'^ - 2r^l\ ' 
i2rf/HiyH - 2ir - 2M)rr - ^ 6M-2. + H(/ + l) j * 

Y sin^ I? 

-4(2rA//)i/2^' - '- ^ , ^ V 7^ H' - 2(r - 2M * - 

r{r^ ~ iMr + + 1) - 2r3 + 6Mr(r - Af) + \2M^ \ ^, , „ , 

r2(6M-2r + W(Z + l)) / ^ 

(3r^?^ + r^/^ - ll^r^ - fr^ + 4Zr^ - ISMr^f + 12Mr^l^ + 6Mr^l* - 24Mr^(? - 1) + 36M^rl{l + 1)+ 

72M2(M-r)) ^ ^ 
(6M-2r + rZ(; + l))^r3 



APPENDIX B: SPIN- WEIGHTED SPHERICAL HARMONICS 

We define the operators c(s) and their adjoints c^(s) by 

c(s) = - ((5 + 2(s + l)aj , ct(s) = (<5 - 2saj , 

respectively. The angular operator defined in (^) takes the form A = c(— 2)c^(— 2). Using the commutation relations 

c{s)c'{s) — c^(s " l)c(s — 1) = — s 

for any real number s we can construct the eigenfunctions of As — c{s)c'' (s) from the the standard spherical harmonics 
y'"*, which fulfill 

0(0)0^(0)^'" = -Iaf'™ = ^l{l + 1)^'™. 

The eigenfunctions of As are called the spin-weighted spherical harmonics, and are proportional to 

^ct(s - l)ct(s - 2) • • • ct(0)r'™, s>0, 



Im — 



c, 



^c{s)c{s + 1) ■ ■ ■ c{~l)Y^"\ s<0 



where 



C'fs = -{l-s + l)il-s + 2)---il + s) 



is a normalization constant. For all s, we have 



AsY^"^ ^ -{I - s){l + s + l)Y^^ . 
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